clc;clear;clf;
p2=[];
s2=[];
for j=3:1000
   for i=3:j
    p2(i)=1/(log(i))^(1025/1000);
   end
   sum(p2)
   s2(j)=sum(p2);
end
plot(s2)